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I. INTRODUCTION 


A. BACKGROUND 


Laser Speckle, that granular, often annoving pattern pro- 
Guced upon laser illumination of a diffusely reflecting 
object, came into prominence as a holographer's nightmare. 
Parle continuing in that role, it is recognized that laser 
speckle, used in double-exposure form, is extremely useful 
for speckle photography. Laser speckle is perhaps the 
easiest form of laser metrology to exploit because in its 
Simplest form, the equinoment required for speckle image 
forming and processing requires only single beam, 
Gouble-exposure, which eliminates the need for vibrational 
isolation. Techniques for reduction of the specklegram have 
been developed successfully and applied to the measurement 
of surface roughness /Il1, 18, 21/7, vibration /6, 20, 39, 51, 
537 and deformation /12, 53/7, and displacement / 4, 7, 19, 


Bete, 537 and strain / 3, 4, 9, 15, 31-34, 40, 537. A most 
convenient form of data reduction for speckle photographs is 
based upon Thomas Young's demonstration in the early 1800's 
of the interference of light emerging from a dual pinhole 


arrangement / 137. Thus a new field of speckle metrology 


evolved. 


Beeecneols OBJECTIVE 
The primary objectives of this thesis were to develop 


various techniques for measuring in-plane and out-of-plane 
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displacement and analysis of strain and vibrations using 
speckle metrology. The objective of using speckle techniaues 
instead of holography procedures in Naval Engineering is to 
enhance the applicability of photographic metrology to the 
industrial environment / 16, 21, 44, 557. Holography, as a 
form of engineering stress, strain, and vibrational metroloay, 
is excellent but requires highly sophisticated equipment, 
close tolerance set up requirements, and a vibrationless 
laboratory atmosphere for exact measurement / 1, 8, 21, 44, 
Be With the development of speckle photography, a single 
beam, easy to setup, high tolerance system of engineering 
metrology is capable of being achieved in an industrial en- 


vironment. / 21, 44/7 
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ii. TREGRY 


A. INTRODUCTION TO SPECKLE METROLOGY 

speckle metrology spans a rather wide variety of tech- 
niques, all of which generally fall into one of two categories: 
direct laser (speckle) photography or speckle interferometry. 
/ 2, 21, 447 The distinction between the two is somewhat 
sudtle since both involve photography and both involve inter- 
ference. That is, speckle photography depends on the optical 
properties of the speckle which involves a condition of inter- 
ference. After vrocessing the speckle photograph, the 
analysis of the interference yields fringes. These fringes 
are directly correlated to the object motion. Therefore, 
there is no great difference between speckle photography and 
speckle interferometry. / 21, 447 Hence, if there are regions 
in the two images where speckle patterns of each are well 
correlated, let the vorocess be referred to as speckle 
photography. If, instead, the fringes form as a result of 
fluctuations in the correlation of the speckle patterns he- 
tween two images, whether or not there is translation be- 
tween the correlated portions of the pattern, let the process 


be called speckle interferometry. 


B. ELECTROMAGNETIC WAVE THEORY 
Light, being a form of electromagnetic radiation, is 
characterized by its amplitude, wave length (or frequency), 


phase, polarization, speed of propagation, and direction of 
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oropagation. (aes, We a3, 537 Light can be described by 
Specifying the temporal and spatial dependence of its electric 
meld vector, E. The simplest tyve of electromagnetic wave 

is the linearly polarized plane wave which has been assumed. 
If such a wave is polarized in the y direction and vrovagates 


in the z direction, the three components of E are: 


E. = 0, 

x 1 
E, = Acos({wt —kz), (1) 
E, = 0. 


where A is the amplitude of the wave. The circular frequency 


w and wave number K are given by 


w= lay 
(2) 
27 


where v is the temnoral frequency and i is the wave 
length. The frequency of light is on the order of 1o-?HzH 
q 
More generally, the direction light travels is described 
> 
by its provagation vector K which has a magnitude K = 2m/A 
and points in the direction of provagation. A plane wave is 
a wave whose phase at any instant of time is constant at all 
; > > 7 ae oe 

Memes normal to K., If r= xi + y} + xk 1s the position 


vector of any point in space, as shown in Fig. 1, the 


equation of a linearly polarized plane wave is: 


x 


ty 
il 


y 


0, 
Acos(wt—k-r), (4) 
0. 


Ml 


? 

? 

. 
Zz 


ty 





A surface over which phase is constant, in this case the 


planes K*r = CONSTANT, is called a wave front. 





Figure 1. Plane wave of light vorovagating in the direction 


Specified by the propagation vector K. From Vest 
[7337 
Lasers emit narrow beams of monochromatic light with 
almost perfectly plane wave fronts. Most lasers emit 
linearly polarized light. Light from a tydical He-Ne con- 
tinuous wave (CW) laser has a wavelength, \X = 632.8 nm. The 
laser emits a beam of about 5 mm. in diameter. Laser light, 
as was mentioned previously, is highly coherent. 
It was also noted previously that the frequency of light 
was on the order of 107° Hz. Practical detectors, such as the 
eye and photographic film, are not capable of responding to 


such extremely ravid variations. Rather, they respond to 


irradiance I which is the time-average energy flux of the 
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light wave. Using electromagnetic wave theory, it can be 
Snown that 


I = eve), / 
Ss 5) 


where € is the electrical permittivity of the medium in 
which light travels; v is the sneed of propagation. The 
key point 1s that the irradiance is proportional to the time 
average of the electrical intensity vector squared; the time 


average 1S indicated by <> ., 


See OLr FRACTION THEORY 


1. Derivation of the Fresnel-Kirchhoff Diffraction 
Formula 


According to the HUYGEN-FRESNEL theory /13/7, the Light 
disturbance at a point P arises from the supervosition of 
secondary waves that voroceed from a surface situated between 
this ovoint and the light source. 


Consider a strictly monochromatic wave 


aba, ,2,t) = U(xX,;7,2) a (6) 


In a vacuum, the space-dependent part satisfies the 
time-independent (Homogeneous) wave equation. 


(y7 + 7) U =O (7) 


where K =uf/c. This is known as the Helmholtz equation. 
Let v be the volume bounded by a closed surface S, 


and let P be any point within the surface; (see Fig. 2), 


2g 








Assume that U vossesses continuous first and second order 
@aretal derivatives within and on the surface. If U' is 
any other function which satisfies the same continuity 


requirements as U, then by Green's theorem 
dU ' dU 
[[[ owe U'V?U)dv [[(e = —U sa) as (3) 


where 9()/dn denotes differentiation along the inward normal 
to S. If U' also is a solution to the Helmholtz equation, 
mee. , if 


(V? + &?)U" = 0, 
OD), 


then it follows that the integrand on the left of equation 


(8) vanishes at every point of v and consequently 


au’ oT 
[[(eg-2 5) ds =0 (10) 





Figure 2. Derivation of the Helmholtz-Kirchhoff integral 


theorem: region of integration. From Born anc 


Wolf /~13 7 
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met U' (x,V,z) = Se which is a spherical wave 
converging to the point P(s=0), where s denotes the distance 
from P to the point (x,y,z). This function has a singular- 
ity at s = 0; and since U' was assumed to be continuous and 
differentiable, P must Le excluded from the domain of in- 
tegration. Surrounding P with a small sphere of radius ¢, the 


integration is expected throughout the volume between S and 


the surface S' of this sphere. In place of equation (10) 
etks e'ks OU] 
Jit} Llea(S)-Fa}s—0 an 


whence 


[[e3(2)-S8]*--[LleF 3) 33) 
--{{ fo (uw) 229 one ane! 
E os 
where 2 denotes an element of the solid angle. Since the 
integral over S is independent of €, the integral on the 
right-hand side may be replaced by its Timi Cinceawe Lue. Oe 
E+ 0. The first and third terms in the integral give no 


contribution in the limit, and the total contribution of 


the second term is 


tks ef OU 
U(P) = ml ios = (> )- — sx} és. (13) 


This equation is known as the Kirchhoff integral 


theorem. It relates the value of any scalar wave function 
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at any point P inside an arbitrary closed surface to the 
value of the wave function at the surface. 

Note that as K+0, the time-independent wave 
equation (9) reduces to Laplace's equation yu = 0, and 


equation (13) goes into the formula of potential theory. 
l d/l 1oU 
BPI a | 7 an (a) — «anh ioe 


Only monochromatic waves will be analyzed, but a 
general form of Kirchhoff's theorem could be derived which 
applies to waves not necessarily monochromatic. 

Consider a monochromatic wave from a point source 
- propagated through an onening in a plane opaque screen, 
and let P be the point at which the light disturbance is 
to be determined. It is assumed that the linear dimensions 
of the ovening, while large compared to the wave length of 
light, are small compared to the distances of both ie and P 
From the screen. 

Homtind the disturbance at P, Kirchhoft*s integral 
is taken over the surface S formed, as seen in Fig. 3, by: 
(1) the opening A, (2) a portion B of the non-illuminated 
side of the screen, and (3) a portion C of a large sphere of 
radius R, centered at P which, together with A and B, forms 
a Closed surface. 


Kirchhoff's theorem then gives 


a 





mmm eI L LL ea(S)-(G)a}s as 


where S 1s the distance of the element ds from P and 
denotes the differentiation along the inward normal to 
the surface of integration. 

The values of U and gU/an on A, B, and C which 
should be substituted into equation (15) are never exactly 
known. However, it 1s assumed that everywhere on A, except 
in the neighborhood of the rim of the opening, U and dU/don 
will not appreciably differ from the values obtained in the 
absence of the screen and that on B these quantities will 
be equal to zero. 


Then Kirchhoff set 








eee. °C 
on w/: U0 = VU"), oe (16) 
aU 
on @: U =0 ee 
where 
Aetkr gui) det ik — =| ee 
Ut) = = - a 1k 5 cos (n, r) (17) 


are the values of U relating to the incident field (see 
Fig. 3b) and A is a constant. The approximations of equa~ 


tion (16) are called Kirchhoff's boundary conditions and 
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meemene Dasis of Kirchhoff's diffraction theory. 





rn ut Pp 
4 
B 
: (b) 
(a) 
fa) Surface of integration (6b) Definition of angles 


Figure 3. Illustrating the derivation of the Fresnel- 
Kirchhoff diffraction formula. From Born and 
Wotton, 


It remains to consider the contribution from the 
spherical portion C. If the radius R is sufficiently large 
so that at the time when the disturbance at P is considered 
no contributions from C could have reacned P, then the 
integral over C will vanish. Neglecting the normal terms 
l/r and 1/s in comparison to k, the Fresnel-Kirchhoff 


Giffraction formula is obtained 


tA eik(r+s) 
O(P) = — yi {[—= [cos (n, 7) — cos (n, s)}dS, (18) 


Closer consideration of the Fresnel-Kirchhoff dif- 
fraction integral, equation (18) , shows, as the 
element ds explores the domain of inteqration, (r+ s) will 


in general change by many wave lengths, so that the factor 
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mec +s), 
a will oscillate rapidly. On the other hand, 


the distances of the points Po and P from the screen are 
large compared to the linear dimensions of the anerture, 
and therefore the factor / cos (n,r) - cos (n,s)7 will 

not vary appreciably over the aperture. All angles are 
assumed to be small. The factor / cos (n,r) - cos (n,s)/ 
may be replaced by coss , where §6is the angle between the 
line POP and the normal to the screen (as seen in Fig. 4). 
Finally, the term 1/rs may be replaced by l/r's', where 

r' and s' are the distances of Py and P from the origin re- 
spectively. The Kirchhoff-Fresnel equation then reduces 


EO 





U(P) ~ — . - { [ew ds. (19) 
Cartesian coordinates are used with the origin i: 

the aperture at 0 and with the x and y axes in the plane 

of the aperture and the positive z direction to point into 

the half space that contains the point P of observation (see 


Fig. 4). 





Figure 4. Diffraction at an aperture in a plane screen. 


Modified from Born and Wolf / d2a7 
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et (Xv 02 ) and x,y,z) are the coordinates of P 
O O 
and P respectively, and ( €,n ) are the coordinates of a 


point Q in the aperture, then 


rt = (9 — 8)? + (Yo — 1)? + 205 one 
gt = (e — 8)? + (y — 7)? + 2, 
rt = 22 + Yor + Zo" (21) 


32 — 72 — y? + 23, 


pez 7!) — (x48 + yon) + 5+ 77, 
g? = 99 — Are + yn) + + 7? (22) 


Since the linear dimensions of the aperture are small con- 
pared to both r' and s', r and s may be expanded as a power 


memeuescetm E/r', n/r', &/r', n/r'. 


rewr’ oS + Yo7 &? + 7? (%yF + Yo)? 


a aaa na as 
ps te Gt? (xf + yn)? 
ewe Oo aa (23) 


Substituting from equation (23) into equation (19) gives 


all 


cos é AI | ikf(é, a Cdn 
a me ve (24) 





a6 





= ae 


where 


{tg eee) ey (ZF + Yo77)? os (25) 


kay Sin, a a 9/3 


ar 23 eC 


The first two directional cosines are (1. 7 To jac (l,m) 


x x 
= — >: Ba 
y 
m= — 2, ey (26) 
and equation (27) may be written in the form 
4 l ee 2 is 2 (27) 
Be) = thy — 1S + (my — mn + ALLS + S| (E+ 97) _ Cot mon _ OS + ma) 
r’ 3’ 


The problem of determining the light disturbance at 
P has been reduced to the evaluation of the integral in 
equation (24). The evaluation can be simplified for two cases 
of particular interest. The quadratic and higher terms in & 
and Nn may be neglected in f£ for Fraunhofer diffraction; when 
the quadratic terms cannot be neglected, Fresnel diffraction 
is being represented. Fortunately, the simpler case of 


Fraunhofer diffraction is of much greater imvortance in 


BeELCS « ro) 
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Vee Prraunnofer Diffraction 

Fraunhofer diffraction occurs when both the incident 
and diffracted waves are effectively plane. This will be the 
case when the distances from the source to the diffracting 
aperture and from the anerture to the receiving point are 
both large enough so that the curvatures of the incident and 
diffracted wave can be neglected (see Fid. 5a). 

If either the source or the receiving point is close 
enough to the diffracting aperture so that the curvature of 


the wave is significant, then Fresnel diffraction occurs (see 
meq. Sb). 


A 


- aN I wW\/ Mn 


perture 


(a) Fraunhofer case (b) Fresnel case 

Figure 5. Diffraction by an aperture. From Cowley /17_/ 
The wave is effectively plane over the aperture if 

the second order and higher terms are negligible. Strictly 

speaking, these terms disappear only in the limiting case 

r'+ ©, s'+o , i.e., when both the source and the point 

of observation are at infinity. It is, however, evident 


that the second order terms do not contribute appreciably 


to the integral if 


Sc 





toi 
o 


1 ] l 
(= sti =) ee (to Sone (1g = ieee) 


r’ 


20 (28) 





The conditions under which equation (28) will be 


satisfied exist when 


(2 + 1 mae 

| p> Et Han and || > Et oes, me 
or if 

St 5=0 and Let, mo'l,m* <a (30) 


The conditions in equation (29) give estimates for 


r' and s' for which the Fraunhofer representation may be 


used. 

In the case of Fraunhofer diffraction, the four quan- 
tities Loe moe Ly and m enter equation (27) only in the 
combinations 


p=l—ly, qg=m— Mp, (31) 


Tne governing Fraunhofer diffraction integral may 


tnen be written 


5° 





Cee of | e~kaE+ an) didn, 
sf 


1 

> EF ={Cl*D 
l /£ 

ear) 


meere c is the constant term in front of the integral of 
equation (24). 

Fraunhofer diffraction patterns for apertures of 
various forms (rectangle and circular, specifically) will 
be studied in detail. 

a. Rectangular aperture 

Consider first a rectangular aperture of sides 
2a and 2b. With origin 0 at the center of the rectangle 


and with n and € axes parallel to the sides (see Fig. 6). 


2a-——_—_——+| 


Figure 6. Pectangular aperture. From Born and Wolf / 13 _/ 


The Fraunhofer diffraction inteqral then becomes 


& 6 5 7 
U(P) = of i e—K(rs+on) d Edn a of” eikps as | e~ikyn dy. 
—a J~b —@ —o ( 3 3 ) 
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Now 


a 
l : 
—tkp§ te eee tna sink 
[ie oS ikp [emia aaa (34) 
with a similar form for the integral of equation (33) taken 
with respect to dn from -b to b. 


Hence the intensity is given by 





sin kpa\? /sin kqb\* (35) 
ur) = (ort (a) (Figs) 
Pac nuinction y = (sin x/x)° is displayed in 
mig. 7, It has a principal maximum of y = 1 at x = 0 
and z@ro minima at x = tn, +27, t37... The roots asymptoti- 
ably approach the values x = (2m + 1)7/2, m being an 
Y 


integer (see Table l). 





Figure 7. Fraunhofer diffraction at a rectangular aperture. 


y = (sin */x)° 
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Within each rectangle formed by nairs of 
consecutive dark lines of the Fourier transform of the 
aperture the intensity rises to a maximum; all their maxi- 
mums are, however, only a small fraction of the central 
maximum of the Fourier transform and decrease rapidly with 
increasing distance from the center (see Fig. 8). The 
larger the rectangular ovening, the smaller the effective 


size of the diffraction pattern of the Fourier transform. 


TABLE 1 
The First Five Maxima Of The Function y = (sin aoe 
2 

x y = (sin x/x) 

0 i 
1.430 = 4.493 0.04718 
2.459 = 7.725 0.01694 
3.470 =10.90 0.00834 
4.479 =14.07 0,00503 


SA em Rg ———E—————E———————————— EE ——————— 
ggg gg a a a I I CE TT ELE ALL ELL ELL LL LEAL LL AI LTC EDL LAI ILD EE I 
ee er rg TT A TT 


= 








feegure 8. Fraunhofer diffraction pattern of a rectangular 


averture 
b. Circular Aperture 
The Fraunhofer diffraction pattern for a cir- 
cular aperture will be investiaqated. Using polar coordinates, 
meee Oo, ©) be the polar coordinates of a particular point in 
the aperture: 


peosd=¢, psind=n; 


(36) 


and let ( w,W ) be the coordinates of a point P in the 
Bmctraction pattern, referred to the geometrical image of 


the source: 


aes = ca) 
weosp=p, wsiny = gq. 


re 
— 
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From the definition of »o and q, it follows that w = ee + a 
is the sine of the angle which the direction (p,q) makes 
with the central direction np =a= 0. The diffraction in- 


tegral of equation (32) now becomes, if a is the radius of 


the circular aperture, 


a (2a 
U(P) = o| { e — ikew cos (0—¥) ndpdé. (38) 
rJ0 


Now using the integral representation of the Bessel 


function J (2): 


i a 
ge cos & et da = JAZ): 


21 0 (39) 


Meauection (38) therefore reduces to 


OP) = 20 | Jub (40) 


Also, the recurrence relationship is 


<. {xnt1 aay = xt] (2), (41) 


giving for n = 0, on integration 


(-aaalede’ = 2h Bis 


0 
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From equations (40) and (42) it follows that 





(43) 
2J, (kaw) 
U(P) = — |, 
oo) cD | kaw 
2 , ; 
where D= 7a. Hence the intensity is given by 
2, (kaw) |? (44) 
as ———— 
KP) = (U(Pye = |", 


where 





I, = C?2D? = ED/i. 


Equation (44) is the formula derived by Airy in somewhat 
different form. / 13 7 

The intensity distribution in the neighborhood 
oz the geometrical image is characterized by the function 
Y = (23, (x) /x)* shown in Fig. 9a. This function y has its 
principle maximum y = 1 at x = 0, and with increasing x, it 
oscillates with gradually diminishing amplitude in a similar 
way to the function (sin ae which was discussed previously. 
The intensity is zero (minimum) for values of x given by 
Jy (x) = 0. 

The minima are no longer equidistant (see Table 2). 
The positions of the secondary maxima are given by the values 


of x that satisfy the equation. 
d = 
ax [oy 00 /x | = 0 
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(a) The function 
205 car 






(b) Airy disk 
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Segure 9. Fraunhofer diffraction pattern of a circular 
aperture. From Fowles / 22 / 





TABLE 2 
The First Few Maxima And Minima Of The Function 
2 
y= (25, (x) 72 








x (23, (x) /x)* 

0 iE Max. 
imezOe= 3.333 0 Min. 
imasse = 5.136 0.0175 Max. 
eee 7.016 0 Min. 
2-679 = 8.417 0.0042 Max. 
ess = 10.174 0 Min. 
3.699 = 11.620 0.0016 Max. 


a 
rnc 





or using the equation Atco) epee lt) (45) 


a ] 
dx i 


by the roots of equations J5 (x) = 0 with increasing 
the separation between two successive minima or maxima 
aoproach the value of wt as in the case of the rectangular 
aperture. 

The results show that the pattern consists of 
a bright disk, centered on the geometrical image at 
Dp = q= 0 of the source and surrounded by concentric 


bright and dark rings (see Figs. 9b and 10). 





Figure 10. Fraunhofer diffraction pattern of a circular 
aperture (the Airy pattern) 6 mm in diameter, 
maanification 50X, mercury yellow liqht i = 


S790A. 
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The intensity of the bright rinas decreases 
rapidly with increasing radius; only the first or second 
rings are bright enough to be visible to the naked eye. 

The bright central area is known as the Airy disk (see Fida. 
10). The Airy disk extends to the first dark ring whose 
size is given by the first zero of the Bessel function, 
namely x = 3.833. The angular radius of the first dark 


ring is thus given by 


IZON L 


Slin§ = a 


(46) 


which is valid for small angles, 6. Here a is the radius 
of the aperture. 

The image of a distant point source formed at the 
focal plane of a camera line is actually a Fraunhofer dif- 
fraction pattern for which the aperture of the lens is the 
opening. Thus the image of a composite source is a suner- 
position of many Airy disks. The resolution of detail in the 
image, therefore, depends on the size of the individual Airy 
disks. If Dis the diameter of the lens, then the angular 
radius of the Airy disk is approximately 1.22\/p. This is 
also the approximate minimum angular separation between two 
equal point sources such that they are just barely able to 
be resolved, because at this angular sevaration, the central 
maximum of the image of one source falls on the first minimum 
Biwtne other (see Fig. 11). This condition for optical reso- 


lution is the Rayleigh criterion. The Rayleigh criterion 
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does not always give an absolute minimum separation for 
resolution, but it is almost always used because of its 


Meeticity. /13, 17, 22, Bey 


Saddle point 





Figure ll. The Rayleigh criterion. 


D. FOURIER TRANSFORMS 
iw CconVvOLUtiONS 

Most considerations of diffraction involve the use 
of the Fourier transform in one form or another. One of the 
most important properties of the Fourier transform involves 
fre concept of a convolution, or convolution integral; 
and for both this and the Fourier transform, it 1s convenient 
to use delta-functions. 

Consider a special delta function, the Dirac delta 
function. A Dirac delta function at x = a is defined by 


QO for x#a 
scr) =| 


ec for x=a 


and (47) 


f d(x-a)dx =1., 


The delta function at x = 0, & (x), can be con- 


Sidered as the limit of a set of real continuous functions, 
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such as Gaussians: 


5(x) = lim | ae exp(-a2s?} | (48) 
As a tends to infinity, the Gaussian function has 

a maximum value tending to infinity and a half-width (l/a) 

tending to zero, but the integral over the function is always 

unity. The delta function can be used as a convenient no- 

tation for any function of integral unity in the form of a 

sharp peak having a width so small it is not experimentally 


meaniticant. 


Mae GQetinition of a delta function is 


(x) = f exp (2aixy}ay , (49) 


imone dimension, the convolution inteoral (or con- 
Moelutions product) of two functions £(x) and q(x) is defined 


as 


Cx) = f(x) #200) =f AX) g(x-X)AX. (50) 


By a simple change of variable. 


———— 


fix) * a(x) =f a(X)fle- XA = a(x) * F(X). (51) 
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The identity operation is the convolution with the 


Dirac delta function: 


f(x) * 5(x) = f(x), 


(5la) 
f(x) * 5 -a) = f(x-a). 


reea ceneral intensity distribution, I(x), each 
Sharp soectral line or each part of a broader spectral line 
will be “spread out" by the "spread-function"” a(x), so that 
the recorded intensity, Tops will be less sharply peaked or 


less well resolved than the original spectrum (see Fig. 12). 


Itx) 
g(x-X) ° 
| | C(x) = I(x) * g(x) 
I(x)-g(x-X) 


Figure 12. A convolution operation. An intensity function 
I(x) representing a spectral line is multiplied 
by the transmission function g(x-X) of a slit, 
centered at X=x; the product of these two functions 
is integrated to give the measured intensity ops (*) > 
Likewise, the blurring of an image due to the im- 
perfections of a camera lens may be described in terms of 
convolution of the ideally perfect image intensity with 


some two dimensional function g(x,y). For a voint source 


a: 





of light, the ideal image would be a two dimensional 


function 46(x,y). The spreading gives 


Lops = 8(x,y) * 5(x,y) = a(x, y) . (rz) 


For a general object, consisting of a large number 
of independently emitting voint sources and having an ideal 


mmage, lL (x,y), 


Ions = G(x,y) «x Ty (x,y) (33) 


Thus each point of the original intensity distri- 
bution is spread into a disk of intensity of light and the 
overlapping of these disks gives a blurring and loss of 
resolution of the image. 

The Fourier transform of a one-dimensional function 


f(x) is defined as 


F [f(x)] = Fu) = J fC) exp {2niux} dx . 
-« (54) 


The inverse transform, pol, is defined so that 


fx) =F FUCOH f Fw) Ro ee (55) 
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Here the convention of including 2™T in the exponent 
is followed. This convention is commonly used in considera- 
tions of diffraction in order to avoid the necessity of 
adding a constant multiplier 27 in equations (54) or Ca 17/ 

2. Diffraction Analysis 

The amplitude of scattering from an object in the 
Fraunhofer diffraction approximation, whether derived from 
the Kirchhoff formulation or from scattering theory, is de- 
scribed by a Fourier transform integral. Thus the diffraction 
amplitude may be described in terms of a distribution in 
Fourier transform space, often referred to as “reciprocal 
space." 

The individual properties of Fourier transforms will 
not be expoundec on here, as they are well documented in 
many references readily available. The remainder of this 
section will demonstrate the application of Fourier trans- 
forms and their relationship to diffraction. 

The amplitude distribution of a very small source or 
the transmission through a very small averture (or slit) in 
one dimension may be described as 6(x) or by &x-a) when it 
is not at the origin. The Fourier transform used to derive © 
the Fraunhofer diffraction pattern is 


F(x) =1, 
ao (56) 


F5(x—a) = exp {2miua} . 


To illustrate the above, the following integral is written: 
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re 








f 6(x—a) exp {27iux} dx . (57) 
Mme integrand is zero except for x = a. Hence the integral 


may be written 


. va (58) 
exp {27iua} yr 5(x~a) dx = exp {2miug} | 
ihe amplitude of the diffraction pattern will be 
proportional to F (u) = F 6(x-a) where u = 1/X . The 


intensity observed will be proportional to F (ia) 2 Le 


ad 


The Fourier transform of a plane wave with respect 


eo t Only, gives 


lee 


g, exp {2mi(v, t—x/A,)} =6(+»,) exp {—2mix/A,}, (59) 


i.e., the delta function in frequency and a plane wave in 


real space. Fourier transform with respect to x gives 


ee 


Fy exp {27i(v, t—x/d,)} = exp {2z7iv, t}8(u- 1A), (60) 


i.e., a Sinuisoidal variation in time and a delta function 
in reciprocal space. Transforming with respect to both 


variables gives 





F, exp {2mi(v, t—x/A,)} = S(vty,)*5(u—3/Ay), 
| (61) 


which is a delta function in both frequency and reciprocal 
space. 

The Fourler transform of all possible waves which 
Can exist in a given medium with respect to both x and t 
is a set of points defining the relationship between the 


frequency v and aS 


(or between the angular frequency 
w and K) which is known as the dispersion relation for that 
type of wave and the particular medium. 


For the translation of an object, the Fourier trans- 


form is written, 


F f(x—a) = FLX) * 8(x—a)] = Fu) exp (2riva} (62) 


Thus the translation of the object in real space has the 
Befect of multiolying the amplitude in reciprocal space by 
a complex exponential. The intensity distribution of the 
Fraunhofer diffraction pattern is given by F(a) = which 
is independent of translation. 
a. Rectangular And Circular Apertures 

The intensity distribution for both the rec- 
tangular and circular apertures is important when studying 
Fraunhofer diffraction. The transmission function of a 


rectangular aperture in two-dimensional form may be written 


U1 
Ui 





Ll if) [x{<a@/2 and [){< 5/2, 
fig) = 


0 (63) 
Then 
a/2 b/2 
F(u,v) = af exp {2miux } dx exp {2mivy } dy af sin(nau) sin(7bv) 
-o/2 —b/2 mau mhy ’ (64) 


so that, for diffraction from a rectangular aperture, the 


intensity distribution is 


> sin?(nau) sin?(mbv) 


‘ 
I(u,v) =a°b° ———- 
a (nau)? (bv) (65) 
; ome 
The maximum intensity at u=vy= 01isab. 
al 


Zeros occur at intervals of a along the u direction, 
which is parallel to the x axis and at b intervals alona 


the vdirection (see Fig. 13). 





F (u,v) 


Figure 13. Diagram suggesting the form of the Fourier 
transform of a rectangular aperture. From 
Cowley / 17_7 
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Thus the intensity function falls off through decreasing 
oscillations in each direction. The dimensions of the dis- 
tribution in reciprocal space are inversely proportional 
to the dimensions in real space. 

The transmission function for a circular 
aperture is written 


1 if (x?+y2)!/? <a/2, 
i (66) 


Ties Ournter transtorm 1S best carried out by 
Writing the Fourier integral in plane polar coordinates. 


The result is: 


J, (mau) (67) 
Pw) -0 





where u is the radial coordinate and J, (x) is the first order 
Bessel function. The function J, (x) Sx is Similar to the form 
(sin x)/x but has a somewhat broader central maximum with 
the first zero at ee ae instead of ane 
b. Single And Double Exposure Specklegrams 

Before calculating the diffraction pattern from 
a double-exvosure specklegram, first consider a single-ex- 
posure photographic transparency of an object illuminated 
by laser light. The recorded image is speckled. If a small 
diameter laser beam is passed through a portion of the 


transparency, the aperture w will be the amplitude trans- 


mission function R(x,y) inside the illuminated area. Thus 


oe 





Us ern x,y) (68) 


and 


rio F {R(x,y)] ne a [Rix,y)] (69) 


Now for a double exposure speckle photographic 
transparency where the object has suffered an in-plane 
displacement Ax in the x direction between exposures, the 


amplituce transmission function is 


Ry (x+y) SO 27) sb oRe x + Ax/2,¥) (70) 
A series expansion of equation (70) yields 


2 te? 
E 25 2 
rere 
where the prime indicates a differential with resnect to 


the variable x. The linear poronerties of a Fourier transform 


provide that 


F Is" (x) | = (iw eee [3 (x) | (725 
where g(x) is the ae derivative of the function g(x). 
Thus 

C=c F[R, (x,y) 

2 4 
iL 

io oc F [R(x,y)]2 2 eh wo + ar ioe + / 

, ‘ (73) 

U=c F iR(x,y)) 2 cos (w, 1 | 





Using equations (72)and (73) and w = K sing 


x 
U = Ws 2 cos ( a K sin? ) (74) 


ret 
l 


I, E tecos (KA x Sind | (75) 


Thus the intensity distribution for a specklegram 
is the diffraction halo from a single-exnosure transvarency 
modulated by Young's fringes. / 3,4/ 

3. Two Dimensional Fourier Transforms 

Two-dimensional Fourier transforms are fundamental 
to all image formation. The Fraunhofer diffraction pattern 
of an object is the Fourier transform of that object. An 
object of low spatial frequency, 1.e., an object that has 
large details, and wide sevaration between details, will give 
a Fourier transform of small elements that are close toaether. 
An object of high spatial frequency will give maxima that are 
wide apart. 

a. Thickness Function 

To find the thickness A(x,y) of a lens at any 
point (x,y) we split the lens into two parts, as shown in 
Fig. 14, and write the total thickness function as the sum 


of the two individual thickness functions, 
ees) = A, (x,y) + A, (x,y) (76) 


Peferring to the geometry shown ier iouce: (| ym tie stale 


ness A, (x,y) is given bv 





Ai(z,y) = An — (Ai - leita y*) 


aay: (77) 
sa = (1-222) 
R,* 


Similarly, 
Ane —(-R. —- VA -— 2? y’) 


rt 7? 
eeu, (1 e Rp . ae) 


where we have factored the positive number R, out of the 


A2(z,y) fae) 


i} 


Beare root. Combining equations (76), (77), and (78), the 


total thickness function is seen to be 


| 2 2 rt yt 
A(ey) = de — R(t - a ned (79) 


A =A A 
where 0 Ol ae 02° 





VR; = x? = y* <x — Ra a VRS — 7 — y? 





Figure 14. Calculation of the thickness FUNCELONe. “som 
Goodman / 25 / 
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b. Paraxial Avyoroximation 
The expression for the thickness function can 
be substantially simplified if attention is restricted to 
portions of the wave front that lie near the lens axis, or 
equivalently, if only paraxial rays are considered. Thus we 
consider only values of x and y sufficiently small to allow 


the following approximations to be accurate: 


2 2 7? 2 








RY? 2R}? 
Fe at r+ y (80) 
i R,? 7 2R:? 


The resulting phase transformation will, of course, represent 
the lens accurately over only a limited area, but this 
limitation iS no more restrictive than the usual paraxial 
approximation of geometric optics. Note that equation (30) 
amounts to approximation of the spherical surfaces of the 
lens by parabolic surfaces. With the helo of these apvroxi- 


mations, the thickness function becomes 


2 ve 1 1 
say) = 4 - >t¥ (2-2) (31) 


c. Phase transformation 


The paraxial approximation to the lens trans- 


formation is: 


= 


2 2 
1 2 


61 





wnere ty (x,y) is the lens transformation function. The 
physical properties of the lens (that isn, Rie and Ro) can 
be combined in a single number f called the focal length; 


which is defined by 
1 1. -) 
fam (x Ra (82) 
Thus the phase transformation may be rewritten 
re 75 . geo aree =) 
 Kz,y) = exp [pkndo] exp | —J 5, y (83) 


This equation will serve as our basic representation of the 
effects of a lens on an incident disturbance. 
One of the remarkable and useful properties of 
a converging lens is its inherent ability to perform two 
dimensional Fourier transformations. 
ad. Fourier Transform Configurations 
Three separate configurations for performing 


the Fourier transform operation are considered (see Fig. 15). 


Objcet 


WI. 


(a) 


(a) Object placed against the lens. 
Figure 15. Fourier Transforming COneL_qumaclens. 
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Object 


| LJ 
a oa 
2) 


(b) Object placed in front of lens. 


Object 
d 
NN 
~~ 
~ 


pe S| 


(c) 
(c) Object placed behind lens. 


mBaegure 15. (Continued). 


In the material that follows, the illumination 
1s assumed monochromatic and the distribution of licht 
amplitude across the back focal plane of the lens is of 
concern. The front and back focal planes of a lens are 
defined as the planes normal to the lens axis situated at 
a distance f in front of and in back of the lens re~ 
spectively. The front of the lens is the side of the lens 
the incoming beam strikes first. 

(1) Object Placed Against The Lens. For the first 
case where the object is placed against the lens, let a plane 


object with amplitude transmittance t, Gy) be olaced imme- 


in 


@eately in front of a converging lens of focal length f, as 
shown in Fig. 15a. The object is assumed to be uniformly il- 


luminated by a normally incident, monochromatic plane wave of 


6 3 





amplitude A, in which case the disturbance incident on the 


mens is 
U,(z,y) = At.(z,y) (34) 


The finite extent of the lens aperture can be accounted for 
Dy associating with the lens a pupil function P(x,y) defined 
by 
PO = 1 inside the lens aperture 
0 otherwise (35) 
Thus the amplitude distribution behind the lens becomes, 


using equation (83), 


Vile) = UdePCea) exp [iE et + v9] (36) 


The constant phase delay associated with the lens transfor- 
mation has been omitted since it does not affect the result 
in any significant way. 

Testing tne distri 5ucron Up (X eve) of field 
amplitude across the back focal plane of the lens, the Fresnel 
@ettraction formula is applied. 
mhus, putting z= f, 

exp i e Gap Se w)| ° 


2f fae 
U,(zy,ys) = pM I] aha ; (37) 


he , 
exp i 5 (x? + “| exp E uf (xz, + wn | dx dy 
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where a constant phase factor has been dropvoed. Substituting 
equation (86) in equation (87), the quadratic vohase factors 


within the integrand are seen to cancel, leaving 


_k 
exp I; of (z,? + w)| 
gf 
Qr 


[f Viz P ey) exp | 3 7 (eer + wn | dz dy 


U;(zy,y7) = 
(38) 


Thus the field distribution Us 1S proportional to the two- 
dimensional Fourier transform of that portion of the incident 
Field subtended by the lens aperture. When the physical 
extent of the object is smaller than the lens aperture, the 


factor P(x,y) may be neglected, yielding 


ok 
A exp 3 af (277+ w)| 


2 i, 
to(z,y) exp | -5 = (an + yi )| 2 d (oo) 
Dio If ie OL ee 


U, (Zy,4s) = 


Evidently the amplitude and phase of the light at coordinates 
(XerVe) are PPeTCereea by the amplitude and phase of the 
Object Fourier component at spatial frequencies 

= = x,/Af f., = y-/*f£). Measurement of the intensitv 
distribution yields knowledge of the vower spectrum of the 


object; the phase distribution is of no consequence in such 


a measurement. Thus 


| [ t.(z,y) exp |-is (zz, + wa | dic dy | (90) 





ieee 
sh0psYs efe 


(2) Object In Front Of Lens. Consider next 


the more general geometry of Fig. 15b. The object, located 
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a distance d, in front of the lens, is illuminated by a 
normally incident plane wave of amplitude A. The amplitude 
transmittance of the object is again represented bv to: In 
addition, let To ery represent the Fourier spectrum of 

the light transmitted by the object, and F a! themrOurler 


ie 
spectrum of the light incident on the lens; that is, 


Fo(fx,fy) = F{Atelan qhilfx.ty) a2 § {U2} conn) 
Filfs.fy) = Fo(fx,fy) exp !—Jdde(fx? + fr)] 
For the moment, the finite extent of the 
lens aperture will be neglected. Thus, letting P = l, 


equation (88) can be rewritten 


exp 3 : ee + y')| 
oF f if 
U;(zy,ys) = eS F; (2 ut) (92) 


Suostitution of equation (91) in equation (92), gives 


exp | =(1- ) 2\ Ge) +] (x u) 


Us (Zi¥s) = pM f \f 


A exp 3 7 x (1 = Far F w)| ff t,(Zo,Yo) CXP | -3 y (tory + vw) | dz,dyo (93) 
baa” = —i >} =a 

Thus the amplitude and phase of the light at coordinates 

(Xer¥e) are again related to the amplitude and phase of the 

object spectrum at frequencies (X- oe, y-/irt ). Note that 

a phase factor again precedes the transform integral, but 

that it vanishes for the very special case do = f, 


(3) Object Placed Behind The Lens. Finally, we 


consider the case of an object that is placed behind the 





lens, as shown in Fig. 15c. The object again has trans- 
mittance toe but is now located a distance d from the focal 
plane of the lens. Let the lens be illuminated by a nor- 
mally incident plane wave of amplitude A. Then incident on 
the object 1s a spherical wave converging toward the back 
focal voint of the lens. 

In the geometrical-optics apvroximation, the 
amplitude of the spherical wave impinging on the object is 
Af/d; the particular region of the object that is illuminated 
is determined by the intersection of the converging cone of 
rayS with the opject plane. If the lens is circular of 
diameter & then a circular region of diameter 2 d/f is illumi- 
nated in the object space. The finite extent of this 
illuminating spot can be represented mathematically by pro- 
jecting the pupil function of the lens alonq the cone of 
rays onto the object vielding an effective voupil function 
P x, (£/4) ry, (£/a)} in the object plane. Note that the object 
transmittance cs will have a finite anerture; the effective 
aperture in the object space is thus determined by the inter- 
section of the true object averture with the projected pupil 
function. If the object averture is fully illuminated, then 
the projected pupil function may be ignored. 

Using a paraxial approximation to the sovheri- 
cal wave that illuminates the object, the field amplitude 


transmitted by the object may be written 
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Af f Jk 
aaa ta) = SP (=. zi Yo ‘) exp |=) od eae + || to(2o,Io) (94) 


Assuming Fresnel diffraction from the object plane to the 
focal plane 
_k 
«A exp E 54 yea: W| 
jad 


ff ew i) & (2. L ue ‘) exp | 7 a (tory + Yoys) | Axo di 
i d d a f Yols o Gio 


Ui (xy/,yy) = : 


(2) 


Thus, up to a quadratic vohase factor, the focal-plane 
amplitude distribution is the Fourier transform of that 
portion of the object subtended by the projected lens aper- 
Eure. 

The result represented by equation (95) is 
essentially the same result obtained previously when the 
object was placed directly against the lens itself. How- 
ever, an extra flexibility has been obtained in the present 
configuration; namely, the scale of the Fourier transform 
is under the control of the experimenter. By increasing d, 
the spatial size of the transform is increased, at least un- 
til the transparency is directlv against the lens (i.e., d= £). 


By decreasing d, the size of the transform is made smaller. /25 7 


BE. INTERFERENCE 
The phenomenon of interference is central to the subject 
Matter presented within. To initiate consideration of inter- 


ference, two different light waves, E, and Eos of the same 
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frequency are superimposed. Since E= E, + E 


l the irradiance 


Oat 
will be: 


I = (E2) = ¢E?) + (ER) + (EE). (96) 


where <> symbolizes the time average or mean. ‘The pro- 
portionality constant ey of equation (5) will be neglected 
for the remainder of the discussion. From the previously 
stated assumption, both waves are linearly polarized in the 


same direction. Then 


E\ = A, cos (wt —ky-r) 


— 


om) 


ore 
E, = A,cos(wt—k,r+¢), ae 


where 6 is a constant relative phase between the two waves, 


and Ky and K, are the wave propagation vectors for Ey and 


E. respectively. Combining equations (96), (97), and (98) 


2 


and carrying out the averaging, then 


= ie i 2 VII; cos 6 


ee) 
where where J, = Aj, 1, = A}, 
(100) 
6 =kyr—ker—¢ 
Bes 2 {Fr — 9 
(101) 
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where @i1s the phase difference between the two waves at any 
location. The irradiance varies from a minimum value 
L 
= — 2 : _ 
Lon I, + L., 2(I,I,) at points where 6 = (2N + 1) 


to a maximum value Tax a elt et 


*5 
D + 2(I,I,) at vnoints 
wnere §6 = 2Nn , where N is an integer. The irradiance 
pattern will be alternating light and dark fringes. 


1. Interference From Two Point Sources Emitting Spherical 
Waves 


A specific important example of interference is that 
of two spherical waves emanating from two point sources of 
light (see Fig. 16). Assume that light is radiated from 
points b and c in phase, that is § = 0. The irradiance at 
any point Ee in space is given by equation (99) with the phase 


difference being f. 
| , P 
fA 
/f 
| IL (\s. - 
7 Zi 


Figure 16. Interference of light emitted by two point 
sources, S, and S.. Modified from Vest / 53 _/ 


If 2 is the perpendicular distance from the axis formed by 


S185 to the point Be and D = d/2 is the distance S1S,, then 
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ry = 2° + {(d - ae (102) 


r5 =\2- ee Giuet x)? (103) 


The phase at oe due to Si and S5 respectively is 
Kr, = 2mxr4/x (104) 
Kr, = 2mro/h (105) 


where kK 2m/xr 


The phase difference generated at the point Pa by 


Sy and S5 1s then 


§ = K(ry = ro) (106) 


Substituting the values for ry and ro 


6 = «Ve cea x)? = V2" aa = 0? | eho 7) 


Assuming the following distance relationships 


ayo, <<) 


perce) 1 (108) 


Factoring 2 out of equation (1907) gives 


2 2 2 2 
(6! 2xd x d 2x X | 
46 K Ll 5 + ; fe 5 lL + 5 5 


Qu 


NO 


using the power identity of 


val 





n n(n-1) x? 
(1 + x) Zo en ee ee (109) 


and neglecting higher order terms 
6 = 2Kxd/ (119) 


mien the locus of points forming a surface of maximum 


irradiance is determined bv setting =e 
where 
meer, = —N 5 ): N= 0, 1, 2 (111) 
1 Z K oe eat oe 


This 1s the equation of a family of hyperboloids 
of revolution about the axis 5,85 connecting the two pvoint 
sources. 

Two light waves which are capable of interfering 
are said to be coherent. Because of the coherence recuire- 
ments, most interference experiments are conducted using 
two images of the same physical source. 

2. Young's Fringes 

A simple division of wave front interferometer is 
used to produce Young's fringes. This is shown in Fig. 17. 
It is simply an opaque screen in which two holes separated 
by a distance D have been cut. This screen is illuminated by 
a point source located a distance Rh from the screen anda 
distance We above the center line axis of symmetry of the 


Screen. The light diffracted by the two holes forms an 


interference pattern which can be observed on a screen placed 


a 





some distance d away. In practice, et y, and D are much 


emailer than 2 and ?., 
Ss Oo 





Opaque screen Ohsrvation 
with two openings screen 


Figure 17. Interferometer to form Young's fringes. 


The irradiance of each of the light waves at y 


therefore will be nearly equal I, = I, = I, SO equation (99) 
becomes 
I= 21, it (Cosa) 


( Bae?) 


i = 41. Boe Cae 2) 


Here tne phase difference is 6 = KA, where A 
is the difference in the distance the light travels from 
the source S to the observation screen at y on the two 


paths. 


ie 
2 2a Z _ 2 = 
AR = b. a (d = Y.) + b, + (d y) 


1 (arr) 


= 
2 eee 2 2 
I + (d+ y.) | b, + (d+ y) 
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Peevyiding by ae and Re to set the equation up to utilize 
the identity of equation (109) and disregarding second 


order and higher terms gives 


2 2 
iday _) (d+y_) 
Moe eo’ Hie 2 Us j|.j, . i S 
s ern 2 eo 
Ss S 
1 (a-y)* 1 (ary)? 
l+s aes i arr eae 
O ‘O 
ae 
= = ie eee 
At Bs 2 eG Q 2 
Dy by 
Ag = - tr 
s O (114) 


Knowing that 6 = KA gives 


Dy, Dy aoe 
= S aS) 
a: 


Where again K = 27 /X is the wave numoer. The 


irradiance is then equal to 


I = 41, cos* ow (116) 
2 


Substituting for $5 into equation (115) gives 
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cos Me) Meany, cae is (alg?) 
Oo yy = _ 


tH 
il 
n 


Soech represents parallel fringe spacing of AL /D. 
3. Scalar Wave Theory 

Further consideration of the subject of diffraction, 
shows the representation of light waves may be modified. 
Seeing that the assumption of unidirectional polarized light 
has been mace, only one vector component of E need be studied. 
In other words, E can be treated as a scalar instead of a 
vector. E is then replaced by the scalar quantity U(x,y,Z,t), 
which is referred to as the optical disturbance. For example, 


a monochromatic plane wave is denoted by 
Mie ec, t) = a(x,y,z) cos(wt —- Ker) (118) 


Similarly, the expression for any momochromatic wave can be 


written in the form 
Wien ,2,t) = a(x,V,Z) cos fut - 6 (x¥ 2) Ciao? 


where a(x,y,z) is the real amplitude of the light wave, and 
8 (x,y,z) is its phase. Of course, the plane wave 1s a 
special case for which 06 (x,y,z) = K*r is constant. Equation 


(119) also can be written as 


U(x,y,2,t) = Re (U(x,y,2) exD (1Wt) 
where 


(129) 


U(x,y,Z) = a(xX,y,2) exp -: a(x. 2) 


io 





is called the complex amplitude of the light. Because the 
temporal frequency of light is so high, there is no need to 
consider the time dependence ext(-iwt) explicitly. The com- 
plex amplitude U(x,y,z) contains all the information about 
the spatial structure of light waves essential to speckle 
photograpny. 

The comolex amplitude of light emitted from a rough 
Surface received at a detector perpendicular to that surface 
1s given by eguation (120). 

The irradiance of the light at the detector surface 
iS a more important quantity because that is what is actually 


detected. The worobahility density function for the irradi- 








ance, Ps i), 1s 
< 1 -I 
Py (I) = Ss OGD ; ) Ci 
<I> 


The irradiance of the snyeckle pattern there obeys 
negative exponential statistics. A measure of the contrast 
of soeckle is the ratio C= o5/ <I> , where I is the 
standard deviation of the irradiance, and I is its mean 
value. For the distribution given by equation (121), the 
sveckle contract is unity. /26/ 

The most probable irradiance in a speckle pattern 


is seen to be zero, that is black (see Fig. 18). 


Hoe oFECKLE THEORY 
The speckle vhenomenon has long been familiar, but only 


the introduction of lasers has brought a deeper understanding 
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and many new applications to this field. /°21 7 Laser 
speckles appear whenever an optically rough surface is il- 
luminated with highly coherent light. The microscopic, 
random height variation of roughness of tne surface needs to 
be of the order of the wave length of light used to illumi- 
nate the object surface or coarser. This encompasses most 
surfaces except tnose highly polished. The optical wave 
Originating from any moderately distant point consists of a 
series of coherent waves, each arising from a different rough 
element of the surface. The optical path difference of these 
various waves may be several wave lengths. These waves 
scattered from the optically rough surface have not only dif- 
fering ovhases but also differing wave amplitudes. The inter- 
ference of the dephased, but coherent, secondary waves results 
in a granular pattern of intensity called speckle. 

1. Speckle Characteristics 

Consider the factors which quantitatively determine 
the characteristics, size, shape, and intensity form of the 
speckle. 

Figures (18 a -f ) are photographs of a tynical 
speckle pattern produced by scattering light from a rough 
surface. 

The physical explanation of speckle is that each point 
on the surface scatters the incoming laser light in such a 
manner that the scattered light from one point interferes with 
that of another point. The randomness of the speckle pattern 
itself is generated by the surface roughness which generates 


phase differences in the reflected waves upon scattering. 
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(b) Speckle size imaged by lens of f/2.8 


Figure 18. Micrographs of sneckle for various lens averture 
settings, showing the sveckle size dependence on 
imaging aperture. Image magnification equaled 
137. Microscopic magnification 300X. 
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Begure 18. 
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Figure 18. 
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To determine the statistical characteristics of the 
irradiance in the observation plane of the sneckle, it is 
convenient to consider an ensemble average at a sincle detec- 
tion point rather than the whole field of view. It is assumed 
Oat the incident light is monochromatic and linearly volar- 
ized, and tnat its polarization is not changed by the scatter- 
eng. 

The most important statistical characteristic of laser 
speckle is its size. From equation (117), it is seen that any 
two voints on the object surface sevarated by a distance 2 
qive rise to fringes of frequency, f = £/Az. The averace 


fringe frequency will be 


(ly Az) (122) 


rh 
il 
Cafe 
rh 
tt 
Wir 


where L is the length of the illuminated area in the y 
direction and z is the distance from the illuminated object 
to the image plane. Hence the irradiance distribution across 


a "typical speckle" will be 
Mmiyye= 1+ cos 2 tT (Ly/3 a) (123) 


where y is the vertical axis value. 


The width of the speckle can be described as the dis- 
tance between points where I drops off to one-half its maxi- 
mum value. This width is 1.5 (A z/%); then a typical sveckle 


width Bs would be 


-— 1.5 (2) (124) 
= Q 
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TO estimate the speckle size in this case, treat the disk 
enclosed by the punil of the lens as a uniformly illuminated 
diffuse surface. If the diameter of the lens is p and the 
image 1s formed a distance z from the lens, D = 1.5 ax , as 
shown in equation (124), while a more rigorous analysis leads 


to a very similar expression 


Pept. 22 Az (125) 


See Fig. (11) for the Rayleigh criterion analysis of resolution. 
If the imaging system is focused on a relatively 


distant plane, z= £ = focal length of the lens; therefore 


b. = 1.22(—&) (126) 


where f/D is the F number of the lens; typical imagings 
systems vary from about f£/1.4 to £/32. See equation (148) 
for a more exact equation for the speckle size. 

The speckle pattern formed by an imaging scattered 
He Ne laser light ( A = 632.8 nm) produces the correspondina 
speckle size which varies from 1 to 24 micron. See Table 3. 
Figures 18(a - £) shows the correlation between F number and 
speckle size. 

For the exverimentation of this thesis, a diffusely 
reflecting object was viewed through a lens or imaging system, 
as indicated in Fig. 19a. 

2. Speckle Photography 
As was stated in the introduction to sveckle metrology, 


whenever diffused laser light is passed through a finite 
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aperture, a speckle pattern is formed. For the science of 
holography, speckle is considered noise because of its 
Memating effect on fringe clarity. /1, 217 

In Fig. 19d, a double-exposure specklegram is seen 
Sera Magnification of 300X. Notice the speckle pairs. These 
speckle pairs developed due to the correlated motion of the 
actual sveckle pattern during the double-exposure recording. 
Processing of this double-exposure specklegram produces a 
Young's fringe pattern from which in-plane translation and 
rotation, strain, and vibrational analysis can be analyzed. 

There are two ovrominent advantages of the speckle 
process over holography. These advantages are simolicity of 
the optical system required for both recording and vrocessing 
the image and the relative ease of interpreting the obtained 


results. 


y 7a 
A 6; 


— 


a 


Object Image (film) 
plane plane 


(a) 


(a) Photographic plate imaging system 


Figure 19. Two-exposure speckle photography for measure- 
ment of in-plane translation. 
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(b) Speckle transparency processing system 





(c) Fringe pattern formed in the image plane of recording 
system. The object was translated in the x direction 
between exposures. 


Figure 19. (Continued) 
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(a) eteTeo eet uze specklegram magnified 300xX. 
Pagure 19. (Continued) 

In-plane displacement will be measured for both 
translation and rotation of a rigid, macroscovicallv-flat 
object using the experimental setup described by Figure 19a. 
An image of the object surface is formed in the film plane 
by a lens of focal length £ and a diameter D. The object 
distance and the image distance 2; are related by the 


thin lens imagery equation: 


Tie 
Laie (127) 


This equation is approximately valid for real lenses and is 
quite accurate when an image is formed using only rays which 
mre close to the optical axis. 

Imagery systems or lenses are specified in terms of 


their relative aperture or F number. The F number is the ratio 


85 





of the focal length f£ to the diameter of the averture D. If 
the same distant object is imaged by a lens of equal F 
number, the average irradiance in their image planes will be 
maentical. 

For an object illuminated by incoming laser light the 
image formed in the film plane will be modulated by a random 
soveckle pattern having a characteristic size De determined by 
the aperture of the imaging lens, as iS seen in equation (126). 


3. In-nlane Translation Measurements (x and y Trans- 
moeton anc ROtation 


If an object translates horizontally by an amount 
L. the relative phase of the light in each of the various 
rays that contributes to the formation of the speckle pattern 
will be unchanged. The entire speckle pattern will translate 
in the film plane of the imaging lens for small translations 
of the object to insure there is no decorrelation of the 
speckle pattern. The sveckle pattern moves a distance ML ler 
the image plane if the object is translated on amount De where 
M is the magnification of the imaging lens. Similarly the 
speckle pattern in the object plane will translate an amount 
=, if the object is translated an amount Lo 

The amount of speckle motion is independent of the 
angle of illumination for in-plane translation. 

To measure in-plane translation a double exposure 
is made. The film is exnosed once prior to the object motion 
and once after the object motion. If the translation is 


Greater than the speckle size be then the resultant negative 
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of the image plane will contain vairs of identical speckle 
patterns as is seen in Fig. 19d. The speckle pairs will be 
separated by a distance ML. The separation ML of the sveckle 
pair can be measured directly by microscopic examination of 
the film. Alternatively, the translation can be measured by 
the optical system seen in Fig. 19a. For this Fourier trans- 
forming system the beam is first collimated. It is then sent 
through the transforming lens. The specklegram is then 
placed behind tne Fourier transforming lens. The collimated 
beam then passes through the transforming lens and the 
specklegram and the speckle pairs recorded on the specklegram 
form Young's fringes in the back focal plane of the trans- 
forming lens as is seen if Fig. 19c. The bright central spot 
seen in Fig. 19c is formed by undiffracted light transmitted 
by the specklegram. Each pair of the corresponding speckles 
act as a pair of identical sources of coherent light which 
form Young's fringes. Since all speckle pairs are separated 
by the same distance ML, all of the Young's fringes overlay 
to form the pattern shown in Fig. 19c. 

Interpretation of the fringes is straight forward. 
As in a Fourier transform, the fringe orientation is normal 
to the in-plane translation L. The magnitude L of the trans- 
lation can be determined by applying equation (117). The 
svoeckles on the transparency are oo ae by a distance Ga; 


the fringes will have a spacing d = a . The in-olane 
S 


translation of the object is given therefore by 
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rAL£ 
Md - (128) 


where A is the wave length of the laser light used to form 
the fringes, Ee is the focal length of the tranforming lens, 
Mis the magnification of the imaging’ system used to record 
the specklegram, and de is the fringe spacing. 
In-plane rotation is the clockwise or counter- 

clockwise rotation of a rigid body. The recording or 
imaging system for in-plane rotation is identical to that for 
in-plane translation. Another type of optical processing 
system must be utilized. Due to the nature of rotational rigid 
body motion every point on the body will undergo a differing 
amount of translation. Therefore, the specklegram must be 
analyzed in a point-by-point fashion. The optical processing 
system for in-plane rotation is seen in Fig. 22a. Ficure 22b 
shows the different fringe patterns that are obtained from a 
counter-closkwise rotation about the plates center. The fringe 
spacing is again inversly proportional to the point motion while 
the fringes are 90° to the direction of the relative motion 
of the body. 

4, Out-Of-Plane Rotation Measurement (TILT) /52, 53/7 

Tilt is out-of-olane rotation. Tilt will be measured 

using a two-exposure speckle photogranvh. The experimental 
setup is identical to that seen in Fig. 19a utilized for 
in-plane translation with a slight modification. Instead of 


focusing the imaging system on the image, tne focus will be 
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placed in the back focal plane or at infinity. Again, a 


speckle vattern with characteristic size b. eee ey lel 


be formed in this plane. 


Film 
plane 


Typical 
speckie 


—— 





Figure 20. Two-exposure speckle photography for measurement 
Cetilt,.© Light rays contributing to the forma 
tion of one speckle are shown. From Vest / 53 7 
Each speckle 1s formed by all light scattered ina 
Sec rcular direction, aS 1S indicated in Fig. 20. If the 
object is translated in its own plane L = ee + LJ the 
speckle pattern remains stationary; however if the object 
tilts by an angle y , as in Fig. 20, the speckle pattern 
in the back focal mane will translate. Each ray contributing 
to the formation of a single speckle leaves the object surface 
at a different location y, but all travel in nearly the iden- 


tical direction. When the object is tilted by an amount jy, 


there is a systematic change of relative phase among these 
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rays. The phase shift can be written 


il 


do =" (ka. —k,)L, (129) 
where K,is the unit direction vector of the i Lanomaina exon 
and Ko is the unit direction vector of the rays contributing 


to the speckle under consideration. For the configuration of 


aug. 20, L = await and 
k,, = 0, k,, = —-sin@, k,, = —cos8, 
= ; (130) 
and fe yO et (131) 
therefore 
C132) 


Ag = = +cos8.)yy. 


Se 
Fe coo ad 


Equation (132) describes a linear variation of vohase 
with y, so the effect is as if all the rays contributing to 
a single speckle tilted by a small angle, (l + cos. )y¥ : 
The corresponding displacement in the focal plane can be cal- 


culated using equation (133). 





Peagure 21, Spatial filtering light of an off axis plane 
wave focused to an off axis spot by a thin lens. 


From vest / 53/ 
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Figure 21 shows what can be considered to be a plane 
wave indicent from an off-axis object point located a larae 
distance away, that is, i ° . This light will be focused 
to an off-axis spot in the back focal nlane. If the plane 
wave propagates at an angle with respect to the optical axis, 


it will focus to a point a distance d- from the axis, where 
de See Can (33) 


We are concerned primarily with the varaxial case; 


hence @ will be small, and 
cane ce = Sin @ = « ise) 


men the translation of the speckle in the back focal plane, 
da, will be 

qd. aes 2. +) Cos oY eisis) 
where fe is the focal length of the transforming lens. tI£ 
the angular anerture of the lens is not too large, the trans- 
lation of the entire speckle pattern will be given by 
Sedation (135). 

To measure tilt, the film is placed in the back 
focal olane of the lens. If a 35 mm. camera is used, the 
camera is focused at infinity instead of on the object. A 
double -exposure speckle photograph is then taken of the ob- 
feet in the infinite focus position. The object is exposed 
once before tilting the object and once after the object is 


tilted. The transparency formed by developing the film is 
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then illuminated by a converging spherical wave formed by 
a lens of focal length Ene as in Figure 19f. This will form 
a sneckle pattern modulated by cosinusoidal fringes of spacing 


d = x Epa. , sO the tilt of the object is given by 


Mr 


Y= Fi +cosd,)d, (136) 


several extensions and applications of the pre- 
viously discussed concepts of speckle photography exist for 
the analysis of vibration and strain. Theory for each of 
these two areas will be discussed here briefly. 

5. Vibrational Analysis / 6, 20, 30, 51, ey, 

In-plane vibrations can be measured by time-average 
speckle photography. / 6, 207 If the object of Fig. 19a 
vibrates harmonically Ze a rigid body, L (t) = (Lx, + Ly) cos wt, 
Where L = (Lx, * = Ly,*) is the amplitude of the vibration. 

A time-average exposure will record speckle streaks of length 
2ML. The speckles spend most of their vibrational period 
near two ends of the streaks; therefore, if the film is pro- 
cessed using the system of Figure 20b, the fringes observed 
in the back focal plane will be similar to those caused by an 
Object translation of 2L. Specifically, the form of the 
fringes that modulate the speckle pattern in the back focal 


plane is described by the square of the zero order Bessel 


Function 





2 ML , 
2a|( Wr \s (137) 
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meee x iS the coordinate parallel to the direction of 
motion. The amplitide of the vibration is determined bv 
measuring the distance d between the dark fringes to the 
feft and those to the right of the bright spot. The fringes 


occur where the argument of J, equals + 2.4, so 


- Mr 
a are (138) 


Ma, 


where A is the wave length of the laser light used to form 
the fringes, =r is the focal length of the fringe imaging 
lens (transforming lens), and M is the magnification of the 
imaging system used to record the speckle photograph. 
Seeserein Analysis /3, 4, 9, 15, 31-34, 40, 53/7 

If the object under study is now strained during the 
experiment, displacement varies from point to voint on the 
object surface. Hence a different method of analyzing the 
speckle photographs must be devised. 

Using the system of Fig. 22, the speckle photograph 
Can be examined one region at a time by illuminating it with 
a small beam such as an unexvanded He-Ne laser bean as shown 
in Fig. 22a. The speckle pattern on the transparency will 


diffract this beam into a diverging cone of angular extent 


aax(L), (139) 


where £/D is the F number of the lens used to record the 


speckle photograph. 


9 3 





Laser 
beam 
a 


Transparency iewi 
(speckle photograph) 





(b) Appearance of fringes at various locations on the 
specklegram of a rigid plate that was rotated about 
a normal axis through the center 0. 
Pigure 22. Formation of fringes by illuminating a double 
exposure photograph with a thin laser beam. 
From Vest / 53 / 
If the small illuminated region contains vairs of 
identical sneckles displaced by a distance da, the Argh c ean 
the resulting pairs of diffraction cones will interfere to 
form a diffraction pattern modulated by Young's fringes. 1B, a 
537 If the fringes are viewed from the photograph at a dis- 
; : AZ 
tance Z from the transparency, their spacing will be d;, = 7 
| 7 S 
Since d. = ML, where M is the lateral magnification of the 
Optical imaging system used to record the speckle vhotograph 


and L is the unknown displacement, 


at Azy 
~ Md,’ (140) 
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Equation (140) gives the displacement of 2OImeEsS: in 
a single small region of the Object surface. The fringes 
are normal to the local in-plane motion. 

An alternative method of analyzing strain was dis- 
covered by Professor Y. Y. Hung. It is the Speckle Shearing 
Camera / 32_/. This method utilizes an image shearing inter- 


ferometer camera, 
CORFRENT 


xX LIGHS SOURCE 
|. (7 AC 
Y Z 


IMAGE SHEARING 
CAMERA 


(a) Schematic diagram of optical setup of an image 
Shearing camera, 
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(b) Imaging detail of the image shearing camera. 


Figure 23. Two exposure speckle photography using image 
shearing camera for strain and vibrational 
measurements. From Hung / 32 7 

The image shearing camera of Fig. 23a has its lens 


divided into two parts by a thin glass wedge. The lens 


images the object, which is illuminated by a diverging beam, 
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inclined at an angle e. to the viewing direction. Since the 
rays Dassing through the glass wedge are deviated, the two 
images, (one focused by each half the lens) are laterally 
sheared with respect to each other as is seen in Fig. 23b. 

Let the plane containing the prism angle be parallel 
to the x-z plane. Rays passing through the wedge are laterally 
shifted in the x-direction. Assuming tnat the wedge is placed 
very close to the lens, the amount of shift dx' in the image 


plane is 


bx’ = 1), (4 = De, 
(142) 


Where D; is the image distance from the lens, UH is the re- 
fractive index of the prism, and a is the prism angle. The 


equivalent shift 4x on the object is 


bx = (ix')/M, (142) 
Where M is the magnification, which is the ratio of the 
image distance D; to the object distance D, from the lens. 


Tnus, 


(143) 


Sy. = Sy (D,/D;) = Da ly - 1) a 


The image shearing camera brings the rays from point 
P(x,y) on the object to meet with that from the neighboring 


point P(x + dx,y) in the image plane. Since the object is 
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coherently illuminated, the rays from the two points inter- 
fere with each other in the film plane. 

The image shearing camera focuses two laterally- 
sheared coherent images in the film plane which interfere 
mith each other. The direction of shearing depends on the 
orientation of the plane containing the wedge angle. Thus, 
by rotating the wedge on the wedge-lens assembly about the 
axis parallel to the viewing direction, tne direction of 
shearing can be varied. The amount of shearing can be 
adjusted, as indicated in equation (143), by varying the wedge 
angle a@ as well as the object distance Do: When the object 
is deformed, a relative displacement between two soints pro- 


duces a relative optical phase change A given by: 


A= 25 (1 + cose ) [w(x48x, y) ~w(x,¥)] » (144) 


sindu(xtdx,y) - ulx,y)| 


where uw and w are the disvlacement components in the x and 
z directions, respectively. If the shear dx is small, the 
relative displacemnts may be approximated by the disnlace- 


ment derivatives and thus equation (144) becomes 
qu 
A ae (1 + cos 8) (ae) + sin BG) (145) 


where the strains are given by 


_ ou _ ov _~ ou 4, ov 
fe. ~ ox Sy ~ ty “xy B ax 
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The term 6x 1s equivalent to the gage length. Should the 
shear be in the y direction, the derivatives in equation (145) 
would be with respect to y. This is achieved by rotating the 
camera lens 90°. tf the object is rotated 90° about the z 
axis, u is replaced by v, the displacement component in the 
7 Glrection. 

By taking a double-exposure recording with the object 
being strained between exposures, a speckle fringe pattern 


depicting A will be generated. Dark fringes occur when 


A = Nt (146) 
where N = De ay, 3S y eee 
TO measure ay this derivative may be isolated 


ox 


using normal illumination (6 = 0). 
However, 1t 1S impossible to isolate du/dx. To 


e ® du e @ 8 ® ‘ 
achieve this as 1solation requires recording two fringe 
du 


patterns using two different illumination angles; then ae 


can be separated from st on a point by point basis using 


the following equation: 


Es 





a ois (147) 


Q> 


X | N, (1 “ cos6,) - N, (1 +. cos8,) | 


sin6, (1 + cos@.,) - sin8, (1 + cos8,) 


where Ny and No are fringe orders corresponding to illumi- 


nation angles 85 and 8. respectively. 
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fee fF Pu THEORY 
1. Film And Optical Requirements 

The sensitivity measurement by sveckle photogranhy 
is limited by the requirement that speckle motion in the film 
plane must exceed the characteristic speckle size, 
be = 1.22 ae -. If £f/2.8 is used as a camera aperture, 
then be = 2.l6ym anda film resolution of about 462 lines/mm 
is required. A high resolution film is necessary for speckle 
photography. $0253 film was used for all speckle experimen- 
tation. This film has a resolution of 2500 lines/mm. The 
maximum frequency of the speckle pattern is FAX = (£)° = 564. 
Since film speed varies inversely with resolution, relatively 
long exposure times were required. 

2. Measurement Limitations 

The upper limit to measurable displacements or tilts 
is set by the decorrelation of the two speckle patterns. MThis 
1s evident because processing of speckle photographs requires 
EOEMatiom of Young's fringes by diffraction of light by nairs 
of identical speckles. A loss of correlation will be caused 
if strains are large enough to alter the microscopic structure 
of the opaque surface under study. Loss of correlation also 
occurs when motions are sufficiently large that the aperture 
of the object forming lens samoles an appreciably different 
portion of the wave front scattered hy the ofject during 
the two-exposures. When this occurs, only a RAG LOmeOr tne 
speckle pattern recorded on the photograph consists of nvairs 


of identical speckles, so the fringes they generate are washed 
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m@aeeoy the light diffracted by the rest of the soeckle. A 
rough estimate of the sensitivity range of speckle photography, 
based on experiments reported, is that motions resulting in 
speckle translations of the order of 15-200um ID een eea el 
plane can be measured. / 6, 20, 537 In the upper limit, the 
object translation equals about ten percent of the diameter 

of the entrance pupil of a typical camera lens at f/4. 

The practical sensitivity range in terms of actual 
object translation or tilt depends on object size, magnification, 
and required depth of focus. / 21, 447 In general, the range 
extends to motions as large as about 1.5 mm. Above this limit, 
direct measure by ordinary double-exposure photography 1s more 
practical. The lower limit of the range 1s about 5um, but 
depth of focus requirements usually raise this limit to about 
100 wm or higher. Hence, in practice there is a large gap 
between the maximum displacement measurable by holographic 
interferometry and the smallest displacement measurable by 


speckle photography. / 44 / 
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Itt. EXPERIMENTAL ANALYSIS AND PROCEDURES 
ee eNO 


fee FILM 

Two types of film were used during this thesis experi- 
mentation. During the Fourier transform, speckle micro- 
photography, and basic Fraunhofer diffraction studies, Polaroid 
mositive negative 4x5 land film with an ASA (55) was used. 

An £/16 aperture and an exposure time of between 1/250 sec. 
and 1/125 sec. were used for optimim results. A Kodak special 
order (S0253) film was used for all speckle photography. A 
Pentax 35 mm camera was used with either a 55 or 85 mm focal 
length lens. For object speckle photography using the 85 mm 
lens with an M = .11 and an F = 2.8 for optimum results an 
exposure time of .5 sec. was used. For object speckle photo- 
graphy using the 55 mm lens with an M = .137 and an F = 4 

and an object distance of 450 mm an exposure time of 2 seconds 
was used for optimum results. For fringe photography only 

the 55 mm lens was used with an F = 2.0 with an object distance 
of five feet giving exnvosure times of 5-7 min. depending on 
the fringe spacing. 

The resolution of the S0253 film is very high (on the 
order of 2500 lines/mm). The higher resolution demanded 
longer exposure times, but this higher resolution was needed 
as will be seen during the discussion of speckle size. 

Lastly, the film developing was done in Kodak D-19 de- 


veloper 70°F for 5% minutes with constant agitation. Then 
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using Kodak ravid fixer, the film was fixed for 3 minutes. 
A film bath followed which usually required 30-40 minutes 
of constant water flow due to the heavy blue tint of tac 


m0253 film. 


Meee oreCKLE SIZE 

As was seen in the theory section under speckle size 
from the governing equation (126) speckle size is a function 
of the wave length of light and the F number of the lmaging 
lens. Table 3 tests the various F numbers and their corre- 
sponding sveckle size, bee when using a He-Ne gas laser as 
was used in the experimentation for this thesis. The wave 


length of a He-Ne gas laser is A = 6328.10 >mm, 








TABLE 3 
Speckle Size Versus Lens F Number And Film Pesolution 


1s ed Z Ze 4 ever 8 ar In 32 


D .x(10- mm) 1.4668 1.544 2.16 3.088 4.323 6.176 8.492 12.352 24.704 


mee On Gen 6647) 462 «323° «231016 «117_~Ss«80 
(lines/mm) 








Looking at the optical system used to record the speckle, 
lens magnification should be included in the speckle calcu- 
lation. / 3, 47 Fig. -24 shows the object-to-image vosition 


which relates the lens equation to the magnification. 
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Figure 24. Formation of a real image with a thin, spherical 
lens. 


Assuming the size of the sveckle dot to be equal to the 


Airy disk in the image plane 


a = 1.22\a si 


9 
: - from equation (125). 


oO 


mere a = 2 in this case. From the lens equation (127) 


~ 


th] 
\ 

It 

wir 


a ab _ 
£ = Tarby or F(a + b) ab 


which relates the speckle size to the image magnification, 


M. It is seen that 


Db. = ae = _ EL ( fi >) 
(148) 
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where 


M = — is the magnification of the imaging system. 

Fig. 25 shows the relationship between the maanification 
M and the sveckle size ba. 

Another area of concern is the effect of maqnification 
change during laser speckle photography. / 457 An assumption 
is made that during the speckle photograph no magnification 
change occurs, where in fact a change does take place. 

For example, in the plane stress case, a magnification 


change arises from the Poisson's effect whereby the plane 


model undergoes lateral thickness changes when loaded in 
its olane. 

In the tilt case, a magnification change occurs due to 
the out-of-plane motion toward the camera thus chandqing the 
object distance and the magnification M. This effect of 
change of magnification is examined in detail by Archbold 
and Ennos / 5, 7, 107 and Dedretti and Chiang / 45/7 The 
problem of magnification change can be minimized by using 
imaging lenses of long focal lengths. 

Mie Airy disk size limitation on speckle size consider- 
ations for resolution is a good approximation, but subpeckle 
size measurements can be effected. / 54 7 This lower measure- 
ment range can easily be reduced to one order of maanitude of 
the usual Airy disk limit. The method of subspeckle size 
speckle photography measurement is discussed in detail by 


Vikram and Vedam /54 7. It requires a known lateral movement 
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being measured. The orientation of the Young's fringes 
gives the unknown deformation. 

There are two other fundamental limitations in the 
recording methods of this thesis. First, laser speckles are 
three dimensional. Speckle patterns decorrelate in a distance 
which is small compared to conventional depth of field for 
an imaging lens system / 3, 4/7. Decorrelation is particularly 
severe with lens apertures wider than about f£/4. Most 
speckle photography of this project was done with a lens 
aperture of f/2.8 therefore decorrelation of speckle was a 
problem. An increased image magnification would reduce the 
problem of the lens aperture effect, but again a low magni- 
fication was used for most experimentation (M = .11). 

Second, speckle photography is significantly vulnerable to 
lens aberrations. Ennos recognized this effect and oublished 
some data with regard to in plane rotation. / 5, 107. A more 
indepth study of the vroblem is given by Stetson / 46 /. 
Stetson's figures showing a comparison of Young's fringes 
achieved by a well-corrected copy lens at unity magnification 
and sveckle halos from a normal camera at the same magnifica- 
tion vividly illustrate the loss of fringe visibility due to 


lens aberations. 


C. FOURIER TRANSFORMS 
In order to better understand diffraction, lens resolution, 
Spatial filtering and other areas of optical physics, a pre- 


liminary study of Fourier transforms was conducted. 
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The first experiment that was conducted consisted of 
comparing differing rectangular aperture spacing. Fourier 
transforms were achieved using Ronchi ruling combinations of 


100-300 lines/mm and 100-450 lines/mm. The results of the 


100-300 lines/mm Ronchi ruling is seen in Ele 26 





Figure 26. Fourier transform of a 100 and 300 lines/mm 
Ronchi ruling combination. 


Three types of filtering were conducted. Using two 
Ronchi rulings of 300 lines/mm, a rectangular aperture grid 
was setup. First, the center of the Fourier transform of the 
rectangular aperture was filtered, which eliminated all lines 
from the reproduced grid pattern as seen in Fig. 27a. 


second , the horizontal lines were filtered out of the 
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Fourier transform of the rectangular aperture, which elimi- 
nated all horizontal lines from the reproduced grid vattern., 
Last, the vertical lines were filtered out of the rourrer 
transform of the rectangulary aperture, which eliminated the 


vertical lines from the grid pattern as is shown in Rilke 2 7D 


J 


(a) Image pattern as a result of filtering the center 
of the Fourier transform. 








(>) Image pattern as a result of filtering the 
vertical lines of the Fourier transform. 


Figure 27. Fourier transform filtering using a rectangular 
Grid . 


Figure 28 shows the lab setup for all Fourier transform 
experimental analysis. A 1.5 milliwatt He-Ne gas laser with 
a A = 632.8 nm was used to produce the monochromatic wave. 

A combination 10X microscope objective combined with a 25 ym 
Spatial filter were used to expand the beam. All lenses in 
this Fourier transform arrangement were spherical double con- 
vex lenses with a focal length fn of 177.8 mm. Another lens 
collimated the beam. The grid was then placed in the optical 
system followed by the transforming lens. The grid pattern 


was then reproduced after passing through the beam filter. 
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A 10X microscone lens was used to expand the beam onto a 
mirror wnich then reflected the reproduced grid image on a 
Frosted glass screen. A magnifying eyepiece enabled the 
Bbserver to view the image grid. 

In order to process the double-exposure svecklegram, the 
Fourler transform system of Fia. 28 was utilized. This system 
of fringe analysis was done using either the film prior to 
or after the transforming lens as is seen in Fig. 15 b& ec 
and described in detail in the Fourier transform section of 


mnis thesis. 


D. IN-PLANE TRANSLATION AND ROTATION 

Figure 29 shows the lab setup for in-plane translation 
and rotation. A 1.5 milliwatt He-Ne gas laser was used to 
produce a monochromatic wave. 

A double shutter was used to facilitate the double-ex- 
nosure speckle nhotogranvhing system. A 40x lens with a 25 um 
spatial filter was used to exovoand tne beam. A 5" x 7" 
2024-T4 Aluminum plate was used for the object. A Pentax 
35 mm camera with an 85 mm lens was used to record the image 
at an object distance of 850 mm. Optimum exposure times for 
this setup were % second for an F = 2.8. 

Howein Dlane translation, a micrometric noider which read 
in 1/10,000 in. and 1 turn = 1/40 in. was used while for 
in-plane rotation a micrometer holder which read in .01 deqrees 
and 1 turn = 0.1 degrees was used. 

Reanstation of .903 in. tO .070 “in. were measured. 


A plot of actual translation versus fringe calculated 
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Laboratory arrangement for recording in-plane 


Figure 29. 





translation can be seen in Fig. 30, while Table 4 shows 
calculation of fringe spacing and fringe density for each 
meranslation. 

The average error for the overall translation data was 
Calculated to be one percent. Both Fig. 30 and the calcu- 
lated error show that there was excellent correlation between 
mae actual and calculated translation. Fidq. 31 (a-f) show 
fringe spacing for various translations. As in the Fourier 
transform experimentation done previously, it can be seen 
that for small translations, speckle pair spacing is small, 
the aperture of diffraction on the specklegqram is small and 
the resultant Young's fringe spacings are large, and con- 
versely for large translation motion, the sveckle pair spacing 
1s large on the specklegram, therefore a large diffraction 
aperture resulting in Young's fringes which are closely spaced 
together. 

Figure 33 shows fringes produced from specklegrams of 
various rotations. Seeing that each point in the snecklecram 
mield is rotated a different amount, individual points on the 
specklegram were inspected. 

As seen in Fig. 32-33, fringes showing almost pure trans- 
iaeion in both the x and y direction were observed at the 
four quarters of the specklegram. At the center of the 
specklegram no relative displacement occurs, therefore Young's 
fringes are not visible. At the outer edge between the 
quarter points of the specklegram, diagonal fringes show the 


existence of both x and y components of the rotation. Moving 
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Calculation Gf Fringe Spacing For Various Translations 


em nnnnnnnnacne nn ncncnn cme n ECE sere eens 


os lines area d, L mm) as on ETFOX * TAG 

HTD er area (mm) mm exp. exo 

e. Ol Se 609 199.09 2250 a. 69 2 loi 5 
fe. 31 8 600 ee) - 3286 1249 0479 
18.45 12 600 50.0 7293 19.4 lon 5 
24.37 16 600 S75 ~ 656 22416 0374 
Bo. 06 18 600 Onto 3 79S 2 oe 19 
Bl. 98 20 600 30 Jos B83 OO 
Bo. 28 24 600 oS ~ 986 ier am 7 Uo 
mo. 3/ LS 365 24,33 1 « OURS Boas moO 2 | 
BO. LZ iy 365 21.47 1.148 45.20 ode / 
m0. 32 les, 365 eo 2 be2e3 9G..0 eS 7 
BS. 51 21 865 7. OCU wecseoZ 95.8 fen 4 1 3 
eee. 95 23% 5 365 15. 332 eee oan, 62.48 0855 
B>.05 25 B65 14.6 1.688 66.47 nozis 
70 26 365 03 lev 3 SoZ 0126 
75 Coo 365 12.807 . beg2ay Vow (6) 06 
Lb = \ Z./Md ) = .6328°107?mn 

exp £ fogs 

Z = 3896.2 mm 
M = .ll 

Mean Error = 3.057% 


114 








(a) Fringes generated due to 0.005 inch translation. 





(5b) Fringes generated due to GieCEO imen translation. 


Figure 31. Fringe spacing for various in-plane translations. 
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(d) Fringes generated due to 0.030 igen tans iacion. 


Figure 31 (Continued) 
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(e) Fringes generated due to 0.040 fic croanSlation. 





(f) Fringes generated due to 0.050 Pee hie tamis Lat] OM. 


Figure 31 (Continued) 
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Figure 32. Double exnosure rotation specklegram showina 
location of various fringe observation mBeints. 
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(a) Lack of fringes observed at Point A. 


Figure 33. Young's fringes produced from various points of 
observation on a rotation specklegram as indi- 


cated in Figure 32. 
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(b) Rotation fringes observed at Point B on 
the specklegram. 





(c) Potation fringes observed at Point C on 
the specklegram. 


Bagure 33 (Continued) 
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(d) Rotation fringes observed at Point D on 
the specklegram. 





(e) Rotation fringes observed at Point E on 
the specklegram. 


Figure 33 (Continued) 
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£) Rotation fringes observed at Point F on 
the speckleagram. 





(gq) Rotation fringes observed at Point G on 
the specklegram. 


Figure 33 (Continued) 
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h) Rotation fringes observed at Point H on 
the specklegram. 





i) Rotation fringes observed at Point I on 
the specklegram. 


Figure 33 (Continued) 
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inward from each of these voints toward the center of the 
specklegram, the relative displacement is less, as would 


be expected, the Young's fringe svacing increases. 


aie OUT-OF-PLANE ROTATION ~ (TILT) 

Figure 34 shows the laboratory arrangement of equipment 
for imaging tilt double-exposure specklegrams. 

The laboratory setupv is identical to the one used for 
recording translation and rotation with two differences. 
First, after the beam is expanded, it is collimated prior 
to illuminating the object. This collimation process is 
necessary for recording all out-of-plane motion using this 
type recording system. Second, instead of recording the 
image in the image plane, the camera is focused at infinity, 
placing the image in the back focal plane of the camera. An 
F number of 2 was used at an object distance of 450 mm with 
an exposure time of 3 seconds per exposure for optimim re- 
sults. A tilt range of 1.5" - 15" was tested producing frinses 
as seen in Figures 35 a and Db. 

Calculations were made of fringe spacing, tilt, and tilt 
angle and compared with the actual tilt observed during 
experimentation. Results of these calculations are tabulated 
mi table 5. 

A graphical comvoarison of calculated tilt angle versus 
calculated fringe is observed in Fig. 36. 

Both the graph, the tabular results, and the mean error 


analysis confirm that there is excellent correlation between 
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(a) Fringes observed due to a tilt of 3" (double 
primes denote seconds of arc). 





(b) Fringes observed due to a tilt of 10". 


Figure 35. Young's fringes produced from various angles 
Geout of Olane rotation (tre) 
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TABLE 5 


Tilt Calculations 








y actual (sec) cd. (mm) Cale. (sec) Error 
aS 6 14.414 so 
10 3 Oe ekg Poe 
5 18 4.804 2039 
4 np Bae O17 
3 a SoS PcG 7 
2 a7 DN SES TI UG 
hye, a2 oS 73 


Mean Frror = 7.70% 


er ee es  ——Ee 


Tilt calculation governing equation and narameters 


y=A2/f(1 + cos 6 ,)d, (in degrees) 


q 


where 


emea2e 10° mm 


Z = 3896.2 mm 


Ff = 55 mm 
z Oo 
Q 30 
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the actual tilt ancle observed during experimentation and 


Ene fringe calculated tilt angle. 


F, MEASUREMENT OF DISPLACEMENTS DUE TO VIBRATION 

Figure 37 shows the laboratory setup for the in-plane 
vibrational experimentation that was conducted in conjuction 
with this thesis. The equipment required for vibrational 
analysis was the same as for tilt experimentation with the 
exception of the plate holders, excitor, power transformer 
and oscilliscope. 

The plate holder was modified to enable small amplitude 
in-plane vibrations to be transmitted to the 5 x 7 aluminum 
plate. A piezoelectric excitor linked in series with a 
spring having a large spring constant enabled vibrational 
amplitudes of .1 - 2 mm range to be transmitted. 

The image was recorded with a 35mm Pentax camera usina 
the double-exposure speckle technique described previously. 
An F = 2.8 with an 85mm lens using an exposure time of 2 
seconds per exposure were used to record the image at an 


image distance of 850 mm and a magnification of paula 


G. STRAIN ANALYSIS 

Figure 38 shows the experimental setup for the strain 
analysis conducted in this thesis. With the exception of 
the plate holder and a modification to the imaging camera, 
which will be discussed later, the equipment setup is iden- 


tical to the tilt experimental equipment system discussed 


previously. 
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The 1lmage shearing camera was used to record the double 
exposure specklegram. In the theory section, the image 
Shearing camera was described in general. Both a 55mm and 
85mm lens were used with the Pentax 35mm camera. The shearing 
prism was a 1° glass wedge with an index of refraction of 1.57. 
With the 55mm lens an image distance of 450mm and an F of 
2.0 with an exposure time of 3 seconds per exposure was used, 
while with the 85mm lens an image distance of 850mm and 
an F of 2.8 with an exposure time of 2 seconds per exvosure 
was used. The image Shearing camera gives a double-exposure 
sSpecklegram which is processed using the system seen in Fig. 39. 

Figure 40 shows the fringe pattern obtained from the 
processed image shearing specklegram. This specklegram can 
be processed in white light without a lens if the shearing 1s 
adequate. 

To produce the fringes of Fig. 40, a 6 in. x 8 in. 
rectangular aluminum olate with its surface sprayed flat red 
was used. The plate was securely clamped along its four 
boundaries. With collimated illumination with a He-Ne laser 
a double exposure of the image was photographed. The plate 
was displaced at the center 0.0005 between exposures. 

A no fringe area was seen during experimentation in areas 
of vlastic deformation. During plastic deformation there is 
no motion of that area of the body, therefore there are no 


local fringes developed in any area where plastic deformation 


1s experienced. 
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Figure 39. Double-exposure specklegram optical vorocessing 
system used to analyze strain fringes. Fram 
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Fringe pattern obtained from the double-exposure 
specklegram for a displacement at the center Or 


the plate of 0.0005 inches. 


Figure 40. 
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Compared with other methods of speckle photography and 
interferometry, the image shearing camera was found to pro- 
duce superior fringes. Besides the imoroved quality of 
fringes there are several other areas where the speckle image 
Shearing camera improves over conventional speckle methods or 
holography. 

First, the setup is relatively simple and does not require 
precision alignment of optics. Second, it does not require 
the special vibrational isolation that accompanies holographic 
measurements. Third, the need of coherence in the illuminating 
light source is Minimized. Although a collimated beam was 
Peamazed, it was not required. Fourth, it provides a wider 
and more controllable range of sensitivity. The sensitivity 
can be changed by either changing the illumination angle 65 
or the amount of shearing. Fifth, the strain is measured 
directly by the image shearing camera. 

While this technique works well for small shearing effects, 
like other speckle techniques the image shearing method is 
highly susceptible to decorrelation when the object under- 
goes large, tilt, translation, or strains. This decorrelation 


deteriorates the fringe visibility. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


mM CONCLUSIONS 


The following conclusions have been drawn from the research 
conducted during the course of this project. 

A comparison of holographic interferometry to speckle 
metrology shows that both systems of measurement have favor- 
able points which support further research in each area. 

The simplicity of equipment arrangement for speckle 
metrology requires no precision alignment of optics as is re- 
guired in holographic measurements Vee Speckle metroloacy 
requires no vibrational isolation system. The use of multi- 
node lasers is possible because speckle metrology minimizes 
the coherent length of the light source required for accurate 
measurements. There is a wider range of sensitivity and an 
ease of sensitivity variation in speckle metrology that is 
fee coundc in holographv 34, 257 BEV Speckle metrology 
requires a much lower film resolution. Decorrelation of the 
speckle pattern due to large tilt, translation, or strains 
deteriorates the fringe visibility. Speckle produced fringes 
are highly susceptible to lens aberration. Even at best, 
speckle produced fringes are of poorer quality than holographic 
fringes. Holographic interferometry presents a much broader 
band of measurement VEE]. Three dimensional motion can be 
studied using holography much more readily than can be 
achieved by speckle metrology / 53 7, 

From the previously listed conclusions, it is seen that 
holographic interferometry used in conjunction with speckle 
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metrology can provide accurate, wide rance, niahlv sensitive, 
means of metrology to Naval Enaineering in general and Naval 
shipyards specifically. 

Specific conclusions derived by actual exverimentation 
are as follows: 

Measurements of in-plane translations were measured for 
a range of 0.001 inches to 0.07 inches with a mean error of 
less than four percent. 

Actual measurements of in-plane rotation were not made, 
but rather a comparison of relative fringe spacing from indi- 
vidual point inspection of the specklegram was performed. 
Direction of rotation and relative size of point displacement 
could be determined from tne resultant fringe patterns. 

Measurements of out-of-plane rotation (tilt) were made. 
The range of experimantation of 1.5" to 15" was too small and 
too few points were sampled to effectively conclude that the 
accuracy of measurement was as precise as experimentation in- 
dicates. Also a more effective method of measuring actual 
tilt should be used to insure the actual tilt measurements 
are exact. 

Vibrational analysis of both in-plane amplitude measure- 
ments and modal shape studies were conducted with limited 
success. The fringes produced from in-plane vibrations gave 
good results, but too small an amplitude range and too few 
points were sampled to make any conclusions concerning the 


effectiveness of this system. 
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Strain analysis using the image shearing camera voroduced 
fringes which showed the strain field over the entire plate 
due to a center displacement of approximately 0.0005 inches. 
Object and image distances can be varied to increase the 
range of shear observation; but this was not done. Therefore, 
the effectiveness of the image shearing camera for a wide 
range of measurements can not be stated. Good fringes were 
achieved using this system for a limited range of shear. 

Each area of study, i.e., translation both in-plane and 
out-of-pvlane, vibrations, and strain, is a broad enough field 


sufficient to support thesis investigation. 


B. RECOMMENDATIONS 

The field of speckle metrology is in its embryo stage of 
Gevelopment. New techniques and methods of analvsis and 
measurement are being discovered continually. The only re- 
striction placed on this field is the imagination of the 
experimenter. 

It is recommended that several of the other methods of 
speckle metrology be investigated. Experimentation in two 
and four aperture Fourier filtering analysis of two and 
possibly three dimensions is being conducted with considerable 
success es pe jae le, BB /.< Dual beam speckle interferometry 
for strain comoonent isolation and vibrational analyses is 
achieving results comparable to that of holography / 15, 39, 


Bye, 51, 53/. 
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In the area of in-plane rotation, it is recommended that 
xX and y component isolation be conducted to measure both 
direction and relative amount of translation using fringe 
spacing and orientation. 

The system of experimentation used to measure tilt was 
found to be highly effective and accurate. Two improvements 
could be made in this area. First, a more accurate method 
is needed for measuring the actual out-of-plane the tilt 
which has occurred. This would validate the error analysis 
conducted in this project. Second, an increased range of 
measurement could be achieved by varying the object and 
image distances. 

The surface of vibrational and strain analyses were barely 
scratched during the course of this project. In conjunction 
with these areas of study, it 1s recommended that speckle 
metrolosy be utilized for beam analysis. Deflection, slope, 
and moment measurements could easily he effected using the 
dual beam / 5,6,10,21/7, Fourier filtering / 21,537, or shear 
camera / 31-337 methods discussed previously. Vibrational 
analysis should be conducted, and the results compared with 
those achieved using holography. 

The area of heat transfer is relatively untouched by the 
literature in conjunction with the rise of speckle metrology. 
It is felt that transient, two-dimensional, heat transfer 
problems could be analyzed using speckle metrology methods. 

Analyses of transient and real time phenomena using 


television speckle recordings have been conducted with 


ears 





relative success by both Maddux / 1 7 and German Da Costa 
of Simon Bolivar University / 21 7. This area of study could 
be undertaken for less than $5,000, greatly expanding the 
field of speckle metrology for the Naval Postgraduate School. 
Currently in industry speckle metrology is being utilized 
for crack detection and non-destructive testing / 21 /, 
pressure vessel inspection / 21 /, and stress concentration 
analysis / 21 7. The applicability of speckle metrology to 
Naval Engineering should be highly evident by the remarks 
entertained in these recommendations and conclusions. Again, 
as mentioned previously, this is a relatively new science; 
the potential for growth is constrained by the intensity of 
the analysts' inquisitiveness. Further studies are con- 
sidered a must, and they are assumed to be very enlightening 


and positive contributors to modern science and Naval Engineering. 
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APPENDIX A 
Notes Added During Proof 
After the thesis was approved, the author continued 
research because of extreme interest in the development of 
speckle metrology. The image shearing camera described in 
the Strain Analysis section of this thesis was used to ob- 
tain slopes of modal amplitudes at various frequencies and 


compare these results with those achieved by Fahey /B-1/. 


I. VIBRATION THEORY 
The specimen of Fig. 41 undergoes a steady state 


Sinusoidal oscillation represented by: 


w(x,y,t) = A(x,y) cos (wt + 6) 


where w(x,y,t) normal displacement of a point P(x,y) at 
a time t. 


A(x,y) = modal amplitude function 


Ww Circular frequency 


Q 


arbitrary phase 
and if normal incidence of light is employed, then A 


is a function of time given by: 
4m aw (x e 
A(x,y,t) = > 6x Oe hey th 
Thus I, (x,y) is a time dependent function given by, 


= f. (x',y') + A(x" yy", t)]] 


s9 








where I, is the intensity distribution received by the 
photographic plate. If the film is exposed to I, (x',y',t) 
for a time T which is large compared to the period of 
vibration, the recorded intensity is an integration of 
I,(x',y',t) Over the time T given by: 


I, (x',y') = = Io(x",y',t)dt 


where I, (x',y') is the integrated intensity. The above 


integration yields: 


mec y ) = T [2 cos<= 2Bx' + £,(x',y') - £2 (x,y) 


sy [32 6x HRC] 


This 1S a case where a high frequency term is amplitude 
modulated by a zero order Bessel function. 


Thus Moire fringes will be formed when: 


0 [Te 6x potted) = 0 


Again, the carrier fringes are not visible by the eye. 
Therefore Fourier filtering is required to make the fringes 
visible. See Fig. 39. 

Figure 41 shows the laboratory equipment arrangement 
for imaging the slopes of the modal amplitudes. The equip- 
ment arrangement is identical to the setup of Fig. 38 for 
Strain analysis with one exception. A piezoelectric excitor 
was placed behind the plate and excited the center of the 
plate at various frequencies and amplitudes using the oscil- 
lator and transformer described previously in the vibrational - 


section of this thesis. 
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Calculations were made previously by Fahey / A-1/7 
from Leissa / A~2/7 for the frequency parameters £, for an 


all sides clamped rectangular plate. 











m n tan 

A i 27.01 
2 i 65.50 
i 2 41.72 
2 2 Poo 


The fundamental frequencies calculated by Fahey / A-1/ 


wer © § 
oe 1 = 327.0 Hz 
Wy > = 506.2 Hz 
9,1 = 794.8 Hz 
"2,2 = 968.45 Hz 


Fundamental frequencies were found experimentally to be: 


4 = 258 Hz 
mize = 548 Hz 
Wo = 752 Hz 
Wo 2 = 1047 Hz 


Modal amplitude fringes were obtained for the three highest 
frequency responses, while at 258 Hz fringes were unable to 
be obtained due to the low amplitude of vibration achieved 


by the excitor. (See Fig. 42a). A dark ring is observed at 
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(b) Slope of modal amplitude at 548 Hz. 


Figure 42. Slopes of modal amplitudes produced at various 
resonant frequencies. 
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(c) Slope of modal amplitude at 752 Hz. 
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(d) Slope of modal amplitude at 1047 Hz. 


Figure 42 (Continued) 
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outer plate where the circular first mode fringes would have 
been seen had a larger amplitude been achieved. 
As observed previously, although the fringes produced by 
holography are much clearer no vibrational isolation is 
needed nor is the precision optical setup required for speckle 
photography. A relatively large range of amplitude variation 
was achieved using the image shearing camera technique. 
Further research is definitely required to enable the 
author to make additional conclusions as to the effectiveness 


of this technique. 
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